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COMMUNICATION METHOD AND APPARATUS FOR MULTI-USER DETECTION 

Field of the Invention 

The present invention relates to a method, apparatus and computer program for the 
5 iterative acquisition of signals for multi-user detection and decoding. 

Background 

Multi-user detection and decoding (MUD) techniques allow detection and decoding 
of transmissions by two or more mutually interfering users over an interference channel. 

In general, MUD involves a detection process in which the received composite 

10 signal is resolved into symbol estimates for each user, and a decoding process in which the 
symbols are decoded to recover their data content, using a forward error correction (FEC) 
decoding algorithm. An optimal joint decoder combines these two processes using a 
maximum likelihood decoding technique to minimize the probability of decoder error. This 
technique has a complexity that increases exponentially with the number of users and the 

15 FEC codeword length, and may not be possible with certain FEC coding techniques such 
as Turbo codes. Hence, it is necessary to devise sub-optimal techniques with reduced 
complexity. 

One sub-optimal approach is to separate the detection and decoding processes. Joint 
decisions are made on each symbol in the detection process, and the symbol streams are 
20 then independently decoded using conventional decoding techniques. 

Another approach is iterative decoding, where soft decisions by the detector are 
input to separate decoders for each user, and the soft decisions by the decoders are fed 
back iteratively to the detector. Examples of iterative decoding algorithms are disclosed in 
*An Iterative Multiuser Decoder for Near-Capacity Communications', Moher M, IEEE 
25 Transactions on Communications vol. 46, No. 7, July 1998 and 'Multiuser Decoding for 
Multibeam systems', Moher M, IEEE Transactions on Vehicular Technology, July 2000, 
Volume 49, Number 4, pages 1226-1234. 

Before detection and decoding can take place, the timing, frequency and phase 
offset of each transmission must be acquired. Timing acquisition is particularly important, 
30 because frequency and phase estimation depend on the correct timing being acquired. 
Moreover, if a signal cannot be acquired or is incorrectly acquired, it cannot be decoded. 
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Statement of the Invention 

According to one aspect of the present invention, there is provided a method of 
iteratively acquiring the timings of a plurality of transmissions in a signal received over a 
multiple access interference channel, the method comprising: estimating the relative 
timings of each of the transmissions; separately soft demodulating, decoding and 
remodulating each of the transmissions to generate soft estimates of each of the 
transmissions; for each transmission, cancelling the soft estimates of the other ones of the 
transmissions to generate an updated estimate of that transmission, and estimating the 
relative timings of each of the updated estimates of the transmissions. 

Other aspects of the present invention include a computer program for performing 
the method, and apparatus arranged to perform the method. 

Brief Description of the Drawings 

Specific embodiments of the present invention will now be described with 
15 reference to the accompanying drawings, in which: 

Figure 1 shows the format of a burst to be acquired in an embodiment of the 
invention; 

Figure 2 is a diagram of multiple transmitters transmitting bursts over a 
multiple access channel; 

20 Figure 3 is a schematic diagram of a multi-user detector and decoder with 

, an acquisition function in an embodiment of the present invention; 

Figure 4 is a diagram of a differential detector for use in the acquisition 
function; 

Figure 5 is a diagram of a coherent detector for use in the acquisition 
25 function; 

Figure 6 is a graph illustrating acquisition performance in a first simulation 
of the embodiment; 

Figure 7 is a graph illustrating acquisition performance in a second 
simulation of the embodiment; and 
30 Figure 8 is a chart showing regions in which acquisition can be achieved in 

the embodiment. 
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Detailed Description of the Embodiments 
Burst Format 

Figure 1 shows one example of a format of transmitted bursts to be acquired in an 
5 embodiment of the present invention. The burst B comprises an initial unique word UW1, 
data D, and a final unique word UW2. The unique words are predetermined bit sequences, 
having low auto-correlation, which are known by a receiver and can therefore be used for 
burst acquisition. The presence of the final unique word UW2 is not essential, but use of 
both unique words improves acquisition performance. 
10 The data D comprises a sequence of modulated symbols x[i], as will be described 

in more detail below. 

A preamble or control word (not shown) may be transmitted before the initial 
unique word UW1, and a guard interval may be left between consecutive bursts in the 
same frequency channel. 
15 As one specific example, the bursts B may be MESP5 or MESP20 packets 

complying with the Inmarsat™ MPDS (mobile packet data service) specification, as 
follows: 
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Table 1 - MESP5 Packet Format 



Modulation 


1 0-QAM 


Input bits per burst 


192 


uoaing rate 


q/7 1 

Off { 


Output bits per burst 


448 


Output symbols per burst 


112 


Preamble 


4 


Initial UW (symbols) 


20 


Final UW (symbols) 


20 | 


Total symbols 15 ms slot 


156 


Guard Time (symbols) 


12 


Symbol rate (ksps) 


33.6 


Slot length 


5 ms I 


Table 2 - MESP20 Packet Format 


Modulation 


16-QAM 


Input bits per burst 


1192 


Coding rate 


y 2 


Output bits per burst 


2384 I 


Output symbols per burst 


596 [ 


Preamble 


4 


Initial UW (symbols) 


40 


Final UW (symbols) 


20 


Total symbols /5 ms slot 


660 


Guard Time (symbols) 


12 


Symbol rate (ksps) 


33.6 


Slot length 


20 ms 



5 Transmitter 

Figure 2 shows a plurality K of users outputting respective bit sequences 
bi[i]...bK[i] encoded by encoders Ci...C K to produce coded sequences di[i]...d K [i], which 
are interleaved by respective interleaves Hi^.Hk to generate interleaved sequences 
E[i(di[i])...IlK(dK[i]), which are in turn modulated by modulators Mi...Mk to generate the 
10 respective sequences of modulated symbols Xi[i]...x K [i] at time L The modulated symbols 
are transmitted in bursts such as shown in Figure 1. The data portion D preferably contains 
an integral number of blocks encoded by the encoders Q...Ck, and the encoders are reset 
after each block, so that the encoding of one burst is independent of the content of any 
other burst. 

15 In one specific example, the encoders Q...Ck are Turbo encoders i.e. parallel 

systematic recursive convolutional encoders, one or more but not all of which have an 
interleaver at the input, as described for example in *Near Shannon limit error-correcting 
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coding and decoding: Turbo codes', Berrou, C, Glavieux, A. and Thitimajshima, P, Proc. 
of ICC '93, pp 1064-1070. The modulators Mi...M K may be 16QAM modulators, as used 
for example in the Inmarsat™ MPDS. 

Multiple Access Channel 

5 The modulated symbol sequences Xi[i] . . .XkW are transmitted on a multiple access 

channel MA such that a set of symbol sequences yi[i]...yic[i] are received at a receiver. 
The effect of the multiple access channel MA may be modelled as: 
y[i]=AWx[i]+n[i] (1) 
where y[i] is the complex vector channel output, 
10 A is a normalised correlation matrix representing the cross-correlation between 

symbol sequences, 

W is a diagonal matrix representing the amplitudes of each user and 
n[i] represents the channel noise. 

Receiver Architecture 

15 The received signal yi[i]...yK[i] is detected and decoded by an iterative MUD 

receiver as shown in Figure 3. A multi-user detector DET takes as its input the output of 
the multiple access channel MA and the current soft estimates (initialised to zero at the 
first iteration) of each user's average contribution to the received signal, subject to the 
current probability distributions on the data. The detector DET outputs updated soft 

20 estimates for each user by subtracting the current soft estimates of all the interfering users. 

The soft estimates for the respective users are soft demodulated by soft 
demodulators DEMi...DEMk, which calculate the posterior probabilities of each possible 
symbol of the modulation constellation. For example, with a 16QAM scheme, for each 
input symbol a probability is calculated of that symbol being each of the possible 16 

25 symbols of the constellation. The corresponding soft detected bits are reordered by 
deinterleavers (not shown, for clarity) and input to soft decoders DECi...DECk which 
refine the probabilities of the coded bits derived from the soft demodulators 
DEMi...DEMk by taking into account the knowledge of the EEC code. The bits are 
reordered once again by respective interleavers (not shown) and output to soft modulators 

30 Mi...Mk, which produce conditional expectations of the coded and modulated symbols 
according to the posterior probabilities calculated by the decoders DECi...DECk. These 
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average symbols are input to a model of an estimated multiple access channel EMA which 
updates the channel estimates for each user, on the basis of estimated channel parameters 
derived by an acquisition function ACQ, and feeds these back to the multi-user detector 
DET for the next iteration. 

5 For each iteration of the MUD receiver algorithm, an acquisition function ACQ 

receives the estimates of each user's contribution to the channel from the multi-user 
detector DET and performs an acquisition algorithm, as will be described below, on each 
of the estimates. The time, frequency and phase detected for each of the users is output to 
the multi-user detector DET for use in the next iteration. At the first iteration, the multi- 

10 user detector DET has no knowledge of the users' contributions to the multiple access 
channel, so its outputs are simply equal to its inputs. The detected time, frequency and 
phase are also output to the estimated multiple access channel EMA. 

Each user's contribution is identified by its acquired characteristics, such as timing 
and optionally frequency and phase, and each separate 'arm' of the MUD receiver operates 

15 on the updated soft estimate for a respective user with the acquired characteristics. As there 
are no current soft estimates during the first iteration of the MUD receiver, no cancellation 
is performed and each arm operates on the same received signal, but with the acquired 
characteristics of the respective user. 

Provided that at least one user is acquired at the first iteration of the MUD receiver 

20 algorithm, the contribution of that user is subtracted from the contribution of the weaker 
users in subsequent iterations, thereby improving the likelihood of successful acquisition of 
the weaker users in subsequent iterations. Hence, the iterative acquisition technique is 
particularly suitable for acquiring weak users in the presence of interference from stronger 
users. 

25 If a false acquisition of a user is made by the acquisition function ACQ, the 

following MUD iteration will attribute a low probability to the decoded signal from that 
user, which will lead to a very low weighting of that user's contribution by the multi-user 
detector DET. Hence, a false acquisition will have little effect on the acquisition and 
decoding of other users. In subsequent iterations, that user may be acquired correctly once 

30 the estimates of the other users have improved. 

The MUD iterations are repeated, on the same received signal yl[i]...yK[i], a 
number of times determined by the desired decoding performance and the acceptable 
processing delay. For example, the number of MUD iterations per received signal set may 
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be fixed at a number likely to give the desired performance under most conditions. 
Alternatively, the MUD iterations may be repeated until the desired decoding accuracy is 
achieved for one or more of the users - this may be determined by the probabilities output 
by the soft decoders DEC1...DECK exceeding a predetermined threshold - subject to a 
maximum number of iterations or maximum processing delay. The decoded bits for each 
user are then output by the MUD receiver. 

The receiver architecture may be implemented in software, programmed for 
example into a digital signal processor (DSP) or other hardware or firmware, which may 
form part of a terminal. The functional blocks shown in Figure 3 do not necessarily 
correspond to discrete hardware components. 

Receiver Functional Description 



15 



20 



25 



Multi-user detector 

The multi-user detector DET has the following inputs: 

a sequence of complex K-vectors (for K users) received from 
the multiple access channel EMA 

a set of sequences of complex K-vectors from the estimated 
channel EMA 



y[i] 



M[i] = 



mti]. 



a sequence of complex KxK matrices. The row vector m[i]k 



is a filter to be used for user k. 
The multi-user detector DET has the following outputs: 

£[i] a sequence of complex K-vectors. The component x k [i]is 

the new estimate of symbol i for user k. 
o k estimated residual interference plus noise variance for user k. 

The multi-user detector DET has the following operation: 
For a synchronous channel, 



f 



x k [i] = (m k [i], 



y[i]-XyjEi] 



(2) 



This is the inner product between the filter for the user k and the results of 
cancelling all the estimated components other than user k from the channel output. For an 
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asynchronous channel, each sequence x k [i] is delayed and equalised to compensate for 
symbol timing and frequency offset, as determined for example by the acquisition function 
ACQ. 

The sequence of matrices M[i] represents a time-varying matrix filter. In the 
5 simplest case, this may be a diagonal matrix of complex values which representing gain 
and phase. These are assumed to be slowly varying and are estimated from unique words at 
the beginning and end of a burst. 

The residual plus noise interference for each user is estimated from the completely 
cancelled signal: 



10 




using a variance estimation technique: 




(4) 



where L is the number of symbols in the bursts. 



15 



Soft Demodulator 



20 



Each soft demodulator DEMk has the following inputs: 

x k [i] sequence of complex numbers corresponding to the 

soft symbol estimates for user k 
42 = (Qo,...Qis) An ordered set containing the normalised complex 

constellation values (assuming QAM, otherwise the 

cardinality will be different) 
{iS* : b = 0,1,2,3} A set containing the indices of the constellation 

points in Q that correspond to bit b being equal to 



zero. 



25 




Noise plus residual interference estimates for each 



user 



the following outputs: 
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(Pk&lPk&lpk&lPk&j) A sequence of 4-vectors (in the case of QAM). 

Element p£[i] is the prior probability, between 0 and 
1, that bit b of symbol i for user k is zero 

and calculates the bit probabilities as follows: 

Pk[i]=Z^ k W-Qj) (5) 

jeS b 

where t| is the zero mean, variance a£ Gaussian probability density function. 
Channel Interleaver 

Each channel interleaver and puncturer receives as input the data and parity bits 
from the corresponding coder Ck and interleaves and punctures them to generate sets of 
bits each corresponding to one symbol for modulation. This type of interleaver and 
puncturer is used in the Inmarsat™ IPDS. Each user uses the same interleaving and 
puncturing pattern. 



15 



20 



25 



Soft Decoder 

Each soft decoder DECk has 
N c 

Ilk, CHANNEL 
Ilk, TURBO 

(p° k [ilpL[ilp^[ilp 3 Jij) 



rck[i] 

and the following outputs: 

(p o k[ilpL[ap^[ilP 3 k[i]) , 



rck[i] 



the following inputs: 

The number of decoding iterations to perform for 

each MUD iteration 

The channel interleaver for the user k 

The Turbo interleaver for the user k 

A sequence of 4-vectors (in the case of QAM). 

Element p£[i] is the prior probability, between 0 and 

1, that bit b of symbol i for user k is zero 

A sequence of scalars representing the prior 

probability that bit i for user k is zero 

A sequence of 4-vectors (in the case of QAM). 
Element p J [i] is the posterior probability, between 0 
and 1 , that bit b of symbol i for user k is zero 
A sequence of scalars representing the posterior 
probability that bit i for user k is zero 
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The posterior coded and uncoded bit probabilities are calculated from the prior 
coded and uncoded bit probabilities using an iterative "soft-in/soft-out" Turbo decpder. 
Techniques for iterative Turbo decoding are well-known in the art, for example: 'Iterative 
Decoding of Binary Block and Convolutional Codes', Hagenauer J, IEEE Transactions on 
5 Information Theory, Vol. 42, No. 2, March 1996. 

Soft Modulator 

Each soft modulator Mk has the following inputs: 

Q- (Qo,..*Qi5) An ordered set containing the normalised complex 

constellation values (assuming QAM, otherwise the 
10 cardinality will be different) 

{Bj : j = 0,1,... 15} A set containing the indices of the bits in j which are 

equal to zero. For example, Bo is the set {0, 1, 2, 3} 
while B s is the set { 1, 3 } 

(Pk&lPktlPk&lpfcH) A sequence of 4-vectors (in the case of QAM). 

15 Element p£[i] is the posterior probability, between 0 

and 1, that bit b of symbol i for user k is zero 

and the following output: 

x[i] a sequence of complex scalars representing the new 

soft estimate of symbol i for user j 

20 

The estimates are calculated as the expectation of the signal constellation according 
to the current symbol probabilities: 

x k [i]=Xe/; (6) 

where the symbol probabilities are calculated according to 

25 Pj= YlpiU\\llQL-piUD (7) 

where B . is the complement set (i.e. the indices of the bits that are ones). 



Acquisition Function 

The acquisition function ACQ has the following input: 
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x[i] a sequence of complex scalars. x k [i] represents the 

new soft estimate of symbol i for user k. 
and outputs the current estimated channel parameters (phase, frequency, timing and 
amplitude). The acquisition function estimates the channel parameters for each of the 
5 sequences corresponding to the different users. 

Estimated Multiple Access Channel 

The estimated multiple access channel EMA has the following inputs: 
x[i] a sequence of complex scalars. x k [i] represents the 

new soft estimate of symbol i for user k 
10 tk, fk> 9k The current estimated channel parameters, including 

time, frequency and phase estimates for each user. 

and the following outputs: 

{^l&l— y K &3} a set °f sequences of complex K-vectors. The vector 

sequence y k [i] corresponds to the current estimated 
15 contribution of user k to the channel output. 

The estimated multiple access channel EMA models the effect of the actual 
multiple access channel according to the current estimates of the channel parameters. For 
example, if the channel is symbol synchronous, 



yji]= 



w lk 
W2fc 



W Kk 



Mi] (8) 



20 where w jk is the estimated complex gain from user k to output j. 

Acquisition Algorithms 

Differential Detection 

In one specific embodiment, a differential detection algorithm is used in the 
acquisition function ACQ. As represented in Figure 4, the algorithm takes the initial and 
25 final unique words UW1 and UW2 of the relevant burst B as input and performs a time 
offset estimation ACQ t . The estimated offset x is then provided as input to a frequency 
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estimation stage ACQf which calculates a frequency offset f using only the initial unique 
word UW1. The frequency estimation stage ACQ f will not be described further. 

In this algorithm, the time offset is detected by the use of differential correlation 
between the received burst and the reference value or values of the initial and final unique 
5 words UW1, UW2; this method is possible because of the low auto-correlation of the 
unique words. 

The unique words UW1 and UW2 may be constant for all bursts, or may be 
selected from one of a plurality of possible unique words - this technique has various uses 
which will not be described here. In the latter case, correlation may be performed between 
10 the received burst and each of the possible reference unique words, and the reference 
unique word having the highest correlation peak is determined to correspond to the unique 
word present in the received burst. 

First, the received signal y(t) is scalar multiplied by the conjugate of the reference 
unique word y(t) with a test offset x which is varied between zero and the maximum time 
15 offset in steps of one sample period: 

s(t) = y(t).y(t-T) (9) 

The power Ro of the result is then calculated as follows: 
R 0 =s(t)xs / (t) (10) 

where s'(t) is the transpose of s(t). Next, the differential power Ri is calculated: 
20 R, =s(t-4T s ^ ol +T^ le )xs , (t + T sampIe ) (11) 

where T sam pie is the sample interval of the receiver and T sym boi is the symbol period. 
For example, in a four times oversampling receiver, T symb oi = 4x Tsampie- 

The ratio of Ri to Ro is calculated for each value of x and the value of x which gives 
the highest ratio is taken as the best estimate of the time offset of the received burst. 
25 The following pseudocode describes the differential detection algorithm: 

Start estimation 

Remove guard and control word from received packet 
FOR 1=0; one sample steps; maximum offset 
30 scalar multiply received signal and conjugate of reference 

UW 

Calculate the power of multiplication results 
calculate differential power of multiplication results 
calculate and record the ratio of the two powers and 
35 corresponding time offset 
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Offset the reference uw by one sample 
END FOR 

Find maximum value of recorded power ratios 
Return the time offset corresponding to maximum value 
5 End estimation 

Differential detection has low complexity when detecting packets with large 
frequency offsets (e.g. > 500 Hz). However, a performance gain may be expected by using 
a computationally more complex coherent detection algorithm, for example as described 
10 below. 

Coherent Detection 

In an alternative embodiment, a coherent detection algorithm is used in the 
acquisition function ACQ. As represented in Figure 5, the algorithm uses the initial and 
final unique words UW1 and UW2 of the relevant burst B to perform coherent estimation 
15 of both time and frequency offset ACQt, f . The estimated time offset x and frequency offset 
are then used by a phase estimation stage ACQq, which calculates a phase offset <p using 
only the initial unique word UW1. The phase estimation stage ACQ<p will not be described 
further. 

First, the frequency offset is estimated by extracting two windows of data, with 
20 length corresponding to that of the initial unique word UW1, from respectively the 
beginning and end of the burst. The first window is correlated with the reference initial 
unique word UW1 and the second window with the reference final unique word UW2. The 
correlation of each window with the reference UW is modelled by: 
r(0 = $(*-*)*y(0 

oo 

25 =jy(Jl-T)ya-T)dA (12) 

where y(t) is the received signal, y(t) is the reference UW and x is the time offset 
applied in the current iteration, and which varies from zero to maximum offset in steps of 
one sample period of the receiver. The symbol * * ' denotes correlation. 

The fast Fourier transforms (FFT) of the correlations are then taken and their 
30 magnitudes summed to yield a vector which is peak picked to find its maximum value. The 
FFT is calculated according to: 



WO 2004/114535 



PCT/GB2003/003288 



14 



*<*)=-^Z*<»)« * (13) 

where N is the number of samples, n is the discrete representation of time t, and k is 
the discrete representation of frequency a>. 

The extraction is done in steps of one sample from zero to the maximum time 
offset. The frequency where the highest peak lies is taken as the estimated frequency 
offset. The sample index of the highest peak corresponds to the estimated time offset. 

The following pseudocode describes the implementation of the coherent estimation 
algorithm: 



10 Start of estimation 

Remove guard and control word from received packet 
FOR 1=0; one sample steps; maximum offset 

correlate received signal with reference unique words 
FFT the correlation result 
15 Find and record highest FFT peak and its corresponding 

frequency 

Offset the reference unique words by one sample 
END FOR 

Find the maximum value among the recorded highest peaks 
20 Return the time offset and frequency corresponding to the maximum 

value 
End estimation 

Optimal Multi-user acquisition 

Under high interference conditions, the single-user acquisition approaches 
25 described above (differential or coherent detection) may limit performance, and it may be 
preferable to use a multi-user acquisition algorithm. Such multi-user algorithms are 
generally too complex to evaluate in detail, but an outline of one possible approach is 
given below. 

The following base band linear Gaussian signal model is assumed. K users transmit 
30 the same, known unique word x(t) to the receiver. The relative time offsets for the users are 
expressed as a vector t = (xi, X2 . . . tr). For simplicity, the data portion D is ignored and it 
is assumed that the signal x(t) is zero outside the interval of the unique word. The vector of 
received signals: 

y, (0 = SV(^^)+ n i (0 ( 14 > 

*=1 
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where is the gain from the k* user's transmitter to spot beam i in a multiple spot 
beam satellite system; ntft) are circularly symmetric complex Gaussian noise with 
covariance E[n t (t)n*j (t + r) J = N 0 if both x = 0 and i = j, and zero otherwise. For simplicity, 
frequency and phase offset are assumed to be zero. 
5 If the channel gains are unknown at the receiver, then an optimal maximum 

likelihood detector makes a joint decision on the gains and the delays xi, x 2 . . . tk, according 
to: 

A,t = arg max p(y I A, t) (15) 

Under the Gaussian assumption, choose the channel gains A and delays t that 
10 minimize the Euclidean distance between the received signal y and the hypothesized signal 
f , where: 

Thus the optimal detector operates according to: 

A,t - arg maxX f |^(0 - y&tfdt (17) 

A,* <ssl " 

t 

15 If the channel gains are known, then only the delays are estimated, i.e. the 

maximization is only performed over t. The maximal likelihood detector essentially 
performs maximal ratio combining on the received signals, involving searching all possible 
combinations of user delays; this is prohibitively complex with currently available 
hardware, although may become feasible with advances in processor power. If phase and 

20 frequency offsets are also present, then all possible combinations of delay, frequency and 
phase offsets for all users would have to be searched. 

Hence, whilst it may be possible to perform multi-user acquisition, single-user 
techniques are preferred. 

Acquisition Simulation Results 

25 The effects of iterative acquisition using coherent detection were simulated in a 

two-user scenario. A 'strong' user has a carrier to interference ration C/I = 10 dB and a 
'weak' user is simulated using a range of smaller values of C/I. The thermal noise power 
Es/No is set the same for both users. Random user data D was created and encapsulated 
into bursts, with random timing offsets between 0 and 0.5 ms, but no frequency or phase 
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offsets, since the aim was to test timing acquisition performance only. The bursts were 
combined in a simulated multiple user channel and additive white Gaussian noise (AWGN) 
was added. 

Tables 3 to 5 below show the results for simulation of 10,000 20ms frames, after a 
5 first and a second iteration of the MUD algorithm, with the weak user C/I = 1, 2 and 4 dB 
respectively. 



Table 3 - Weak User C/I = ldB 



Es/No 


Strong User 
Errors 


Weak User Iteration 1 
Errors 


Weak User Iteration 2 
Errors 


0 


1 


3100 


31 


1 


0 


2110 


0 


3 


0 


1500 


0 


6 


0 


990 


0 



Table 4 - Weak User C/I = 2dB 



Es/No 


Strong User 
Errors 


Weak User Iteration 1 
Errors 


Weak User Iteration 2 
Errors 


0 


1 


990 


8 


1 


0 


770 


0 


3 


0 


310 


0 


6 


0 


110 


0 



10 

Table 5 - Weak User C/I = 4dB 



Es/No 


Strong User 
Errors 


Weak User Iteration 1 
Errors 


Weak User Iteration 2 
Errors 


0 


1 


790 


5 


1 


0 


220 


0 


3 


0 


46 


0 


6 


0 


6 


0 



Figure 6 shows the acquisition performance for the weak user before and after 
(shaded region) the strong user was subtracted. These results show that the iterative 
15 acquisition method can successfully acquire users that cannot be acquired by single-user 
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non-iterative techniques. For Es/N 0 > 1 dB, the weak user can be acquired (error rate less 
than 10" 4 ) using the iterative method. 

A further simulation was performed using Inmarsat™ BPDS bearer sub-type L3, 
which has a coding rate of 1/3, under fading conditions (C/M = 15 dB, fading bandwidth 
5 20 Hz, multipath delay = 0). Otherwise, the parameters were the same as in the simulation 
above. The results are shown in Tables 6 to 8 below. 



Table 6 - Weak User C/I = ldB 



Ej/No 


Strong User 
Errors 


Weak User Iteration 1 
Errors 


Weak User Iteration 2 
Errors 


1.25 


1 


4360 


4 


2.25 


0 


960 


0 


3.25 


0 


600 


0 



10 Table 7 - Weak User C/I = 2dB 



Es/No 


Strong User 
Errors 


Weak User Iteration 1 
Errors 


Weak User Iteration 2 
Errors 


1.25 


0 


2810 


2 


2.25 


0 


1300 


0 


3.25 


0 


525 


0 



Table 8 - Weak User C/I = 4dB 



Es/No 


Strong User 
Errors 


Weak User Iteration 1 
Errors 


Weak User Iteration 2 
Errors 


1.25 


1 


1100 


2 


2.25 


0 


620 


0 


3.25 


0 


320 


0 



The results are shown in Figure 1. Acquisition is again possible in regions where it 
15 would not be possible using single-user non-iterative techniques only. For Es/N 0 > 2.25, 
the weak user can be successfully acquired. 

Figure 8 shows the region in which the iterative acquisition method can be used 
successfully to acquire a weak user in the case where the strong user has C/I = 10 dB. The 
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horizontal axis represents the C/I for the weak user, while the vertical axis represents 
signal to noise ratio (Eg/No), which is common to both users. The shaded region above and 
to the right of the line represents the conditions where acquisition of the weak user is 
achievable with error rate less than 10" 4 . Hence, the iterative method greatly improves the 
5 acquisition performance for weaker users. 

In summary, it has been demonstrated that the use of iterative acquisition 
techniques in combination with iterative decoding techniques gives superior performance 
to non-iterative initial acquisition techniques in a multi-user detector. The superior 
performance is due at least in part to the use of information derived by a decoding iteration 
10 to improve the performance of the next acquisition iteration. This advantage is not 
dependent on the specific details of the examples and simulations described above. Hence, 
the skilled person will recognise that the present invention is not limited save as defined in 
the following claims. 



15 



